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We study the short-time and medium-time behavior of the survival probability in the frame of 
the A'^-level Friedrichs model. The time evolution of an arbitrary unstable initial state is determined. 
We show that the survival probability may oscillate significantly during the so-called exponential 
era. This result explains qualitatively the experimental observations of the NaT decay. 



I. INTRODUCTION 



OO i Recent developments in femtosecond laser optics, see for example the XXth Solvay Conference on Chemistry [||, 
opened new possibilities for the study of quantum transitions, which are a very important subject of the quantum 
theory. In a series of works, Zewail et al |^|-^ applied femtosecond transition-time spectroscopy for the probing 
of chemical reactions. Following the work of Kinsey et al. they attempted in paper to track wave packet 
trajectories in the dissociation of Nal. 

The shapes of the ground state potential for Nal and of the quasi-bound potential of the Na(^Si/2)+I(^P3/2) system 
I suggest a mechanism of the induced dissociation process. The femtosecond laser pulse brings the Nal molecule to 
the state of quasi-bound ions. The distance between the ions reaches the region where two potentials have similar 
values due to vibrations of a Nal excited state. Then the transition from Na(^Si/2)+I(^P3/2) quasi-bound state to 
Nal continuum state occurs resulting in the dissociation of the molecule. 
' After an initial exciting laser pulse, the experiment shows oscillations of the Na(^Si/2)+I(^P3/2) population, which 
Oh! are explained in |^ by wave packet propagation. The direction of the wave packet propagation is correlated with 
^Ij • the oscillation (extension and contraction) of the Nal bond. The quantum dynamics calculations are based on a 
C [ time-dependent perturbation formalism. 

This problem is an example of the interaction of the discrete spectrum with the continuous spectrum, which 
was extensively discussed in the literature starting from the work of Friedrichs j?). Indeed, the energy states of 
Na(^Si/2)+I(^P3/2) are the excited state embedded into the continuum states of the decay products. Therefore, the 
time dependence of the Na(^Si/2)+I(^P3/2) population is described by the survival probability of the excited state 
, prepared by the laser pulse. 

The original Friedrichs model ||^ contains two discrete energy level, a ground state and an excited state, coupled 
with the continuum, being bounded from below. The time dependence of the survival probability of the excited state 
has been studied both theoretically |p|-p^ and experimentally . It is exponential with a short non-exponential 

initial era and a non-exponential long tail. As a result, Friedrichs models are very appropriate for the discussion of 
the particle decay and for the description of dressed unstable states ^-po| . The analytical structure of the A^- level 
Friedrichs model has been widely discussed pl|-[2^ , and the possibility of the oscillations of the survival probability 
was pointed out in 

In the present paper we shall show that the A^-level Friedrichs model can also explain the oscillations of the survival 
probability of the excited state observed by Zewail and co-workers . Several excited levels are necessary in order to 
construct a wave packet, which can exhibit localization and nonconventional time evolution. In Section 2 we present 
the model and describe the exact solution diagonalising the Hamiltonian. Using the relation between eigenstates of 
the unperturbed Hamiltonian and the total Hamiltonian, we describe in Section 3 the time evolution of the basis 
states. Specifying the formfactor of the interaction, we show in Section 4 the appearance of oscillations already for 
the two level Friedrichs model. In Section 5 we demonstrate that the survival probability of unstable states in the 
A?^- level Friedrichs model is in fact very close to the one obtained in the experiment jSj. 
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II. MODEL AND EXACT SOLUTION 



The Hamiltonian of the Friedrichs model generalized to A^-level is: 

H = Ho + XV, (1) 

where 

N 



Hq ^^^uJk\k){k\ + / dujuj\u!){uj\, 
k=i { 

V =J2 [ dcjfk{u;)m{u;\ + \u;){k\). (2) 

Here |fc) represent states of the discrete spectrum with the energy ujk, cuk > 0. We assume the simplest case that 
Wfe ^ uJk' for k 7^ k'. The vectors \uj) represent states of the continuous spectrum with the energy w, are the 

formfactors for the transitions between the discrete and the continuous spectrum, and A is the coupling parameter. 
The vacuum energy is chosen to be zero. The states and jw) form a complete orthonormal basis: 

{k\k')^dkk', {lu\lo') ^Siu^-u^'), (c^|fc)=0, k,k' = l...N, (3) 



Y,^k\k){k\+ / dL0\uj){L0\=I, (4) 

fe=l { 

where 6kk' is the Kronecker symbol, 5{uj—uj') is the Dirac's delta function and / is the unity operator. The Hamiltonian 
Hq has the continuous spectrum on the interval [0, cx)) and the discrete spectrum lui, Wfc embedded in the continuous 
spectrum. 

As the interaction XV is switched on, the eigenstates |fc) become resonances of H as in the case of the one-level 
Friedrichs model |0] . Let us consider the eigenvalue problem for the iV-level Friedrichs Hamiltonian (0) 

H\^S^) =iu\^'^). (5) 

We shall look for the solution of Eq.(|^) in the form: 

oo 

1*^) - 5]^feM|fc> + J dcj'i;{iu,u;')\u;'), (6) 
k 

where ipk{i^) and ip{uj,uj') are unknown functions. Inserting into (|^) and making use of the orthogonality relations, 
we obtain for them a system of equations: 

oo 

(wfe - cj)-0fe(w) + A / (itj'/fc(w')7/)(w,tj') = 0, 

(7) 

k=l 

Eliminating tp(uj,uj') from this system, we arrive at the following equation for ipkiuj): 

N 

Y,Gll,{Lo)^k'{io) = -CXh{Lo), (8) 
fc'=i 

where C is an arbitrary constant. Gkk'{^) are the matrix elements of the partial resolvent which is: 

DO 

ri-l ( \ I \x \2 /"j I fk{^')fk'(yj') 

^fcfe'(^) = - ujjSkk' - A duj ; . (9) 
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Under certain conditions (which will be specified below, see also 26|), the function Gkk'{z) is analytic everywhere in 
the first sheet of the Riemann manifold except for the cut [0,oo). In this case, the Hamiltonian H has no discrete 
spectrum. The solution of Eq. (p|) is given by 



N 



M^) = ~C\ J2 Gkk'{io±iO)fk'{u;). 

k' = l 

With this equation we find 'iJj{uj,uj') from the system (|^): 
ipiuj,u;') = C 



(10) 



„ X^j:k.k'=iM^')Gkk'{u;±tO)h'{u;) 
d(UJ — uj ) ^ 



w — w' ± zO 



(11) 



The eigenvalue problem (||) has two sets of solutions 

N f °° 



out 



W) - \k) 



(12) 



which correspond to the "in" and "out" asymptotic conditions. The value C — 1 corresponds to the orthonormalization 
condition 



in in ^5{i^-J). 

out out 

We can also prove the completeness condition 

OO OQ 

/da;|*^>„, „, (*^| - Vwfc|A:)(fc|+ j duo\uj){uj\. 

J outout J 


Hence the new states diagonalize the total Hamiltonian (||) as 

OO 

J outout 


The proof of completeness is based on the matrix formula 

A-^ - B-^ = A-^{B - A)B-^, 

from which we can derive: 

N 

Gkk'(yi + - Gkk'{^ - i^) ^'2.m>? ^ G'fc;(cj + iO)/;(tj)/„i(cj)Gmfe' (tj - iO) , 

1,711 — 1 



Using the asymptotics: 



we prove other useful relations for G: 



Gkk' (t^ )u^=^. : + O 



UJ — UJ' \0J — OJ' 



N 



Gkk'K'^ ± zO) = A I duj fi(uj )fm[i^ ) — 



l,m—l 



and 



N 



duj fi{Lj)fjn{uj)Gki{uj + iO)Gjnk'{^ - iO) ^ Skk' 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



l.m—l 



Due to the completness of the new basis ( p^ ) the old basis vectors may be expressed in terms of the new ones as: 

oo oc 

\k)= J dLo\^^),^,^{^^\k), \lo)= j duj'\^^,),^U'^M (20) 



where \n{^u)\k) and \n{'^uj'\^) are the complex conjugates of {k\^^)in and {'jjY^u)')^ respectively, which may be 
obtained from (p^): 

N 

{k\^J)n. = -A^/;HGfci(^ + «0), (21) 

1=1 

I ^ s( \- \'h{u)h{u')G-,,{u') 

(W|*..')in = 5{UJ - UJ ) ~ uj'^LO~iQ ■ ^ ^ 

Inserting complex conjugate of (|2^) into (|2^) we obtain the inverse relation in the form: 

N °° 

\k) = -A^ I dLofi{uj)Gki{uj - iO)|*^)in (23) 
1=1 i 

^ — ' / UJ' — UJ — l\J 

k,l = l Q 

These inverse relations will be used for the calculation of the time evolution of \k) and \lu) in the next section. 

III. TIME EVOLUTION 

Using the known evolution of the state |^a;)in, 

we can find the evolution of the eigenstates of Hq: 

N °? 

\k)t = -XYl duje''^'Muj)Gki{uj~tO)\^^)in, (25) 
1=1 { 

1^)* = e--*|vl/j,, - XUu) ( du;-e— '^ ^^'l^'^^^'^V ^Oin • (26) 

^ — ^ / i,y — in — 7(1 



Using (|12|), we obtain the representation 

N N 



|fc)t = ^Afci(t)|0+A^ f dLufi{to)gki{to,t)\L^), (27) 

1=1 1=1 g 

N °° 

\uj}t = e \uj) - X y fi(uj) / duj fk{uj) 

/ UJ — UJ 



N 



J2 \fk{oj)9M{uJ,t)\l) (28) 



k,l=l 



in terms of the time-dependent matrix functions Aki{t) and g{ijj,t): 



N 



Akiit) = ^ / duje-''^'U{Lj)fn{u;)Gkm{u; + iO)GUlj - iO)\l), (29) 



gM{uJ,t)^-e-'^'Gu{ou~tO) (30) 

N °? 



AT ^ 

" ,^_,^'t/™(^(c^')Gfc™(c^'-iO)Gz„(c^' + zO) 



oj' — uj + iQ 

With the help of (116^, we can rewrite (|29|) in the form 



oc 

Aki{t) - 2^ / dioe-'"' {Guiio + tO) - Gki{oJ - tO)) = Aoje-'^'Guioj) , (31) 

c 

where the contour C is shown in Fig. 1. With the help of (p^), we rewrite ( pO[ ) in the form 

AT ^ I 

5fe,(c^,i) = A2 V duj'f„,{uj')UJ)Gkm{uj' -iQ)Gin{uj' + — - . (32) 

J LO' — U) + lO 

m,n— 1 Q 

The integrand in ( |3^ ) does not have any singularity at uj' — oj, therefore iO in the denominator becomes redundant. 
Then using ( [T6| ) we obtain 



oc 

iu t ^—iujt 



1 f e"*"^ — e 

9ki {u:,t) = — dcj' {Gki {cj' + ^0) - Gki (co' - lO)) 



2-Ki J LO' ~ UJ 

p — iLJ t ^ — iut 

duj'Gkiii^') ^- , (33) 

c 

where the contour C is shown in Fig. 1. For real lo > Q the term with the factor e~*"* vanishes because it does not 
have any singularities outside the positive part of the real line. Then we have 

1 r -iuj't 
gkiiio,t) = — / dcu'GMiio')- . (34) 

c 

One can easily check the following relation between Aki{t) and gki{uj,t): 

Aki{t) ^ [i^- uj^ gki(.'^,t) . (35) 
The time evolution of a state |$), which is a superposition of the eigenstates of iJo, 

N 



|$)=^afe|fc), (36) 



may be obtained with the help of (p 



k=l 



N 



|ci>(t))=^a,|fc)*. (37) 



fe=i 

The survival amplitude A{t) of this state is 

N N 



A(t) = ($|$(t)> = akal,{k\k')t^ ^kal,Akk'{t). (38) 



k,k' = l k,k' = l 



Changing the contour of the integration C to Ci in Akk'{t) as shown in Fig. 1, we arrive at 
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1 



(39) 



where r;[^., is the residue of Gtk' (w) at the pole 



^kk' = resGfcfc'(w)|, 



(40) 



The first term in ( |39| ) corresponds to the contribution of the poles Zj while the second term is the background 
integral, which gives rise to so-called long tail behavior [§,0. It is known that the integral term plays essential role 
for very long as well as very short times. In the case of very short times we have the well-known Zeno and anti-Zeno 
regions |^,^ 14 1^- If we consider the intermediate "exponential decay" era, the integral term can be neglected 

because in this time scale, it is of the next order in compared with the first term. 

The same result for Akk'{i) (|3^) is obtained in Appendix A in terms of Gamov vectors (A14). In the intermediate 
"exponential" era, the main contribution to the survival probability comes from the Gamov vectors as one may neglect 
the integral term arising from the background. 



IV. TWO LEVEL MODEL 



The rich structure of the model involving more than one level, will be first illustrated with example with two excited 
levels by choosing the formfactor in the form similar to 



Ik{^) = — 



For this formfactor the matrix element Gj,^, {uj) (||) is 



7rA2 



Pk + Pk' (v^ -I- ipk){V^ + ipk') ' 



(41) 



(42) 



where the first sheet of the complex to plane corresponds to the upper half of the complex y/uj plane. The square 
root is defined with the cut [0, -f oo) such that ^/uj > at the upper rim of the cut. For pk > all singularities of the 
integral in expression (0) are on the second sheet. In the case of two levels the matrix is 



/ 



{uji - Cj) + 



7rA2 



ttA^ 



\ {pi + P2){V^ + ■ipi){V^ + ip2) 
The 2x2 matrix representing the partial resolvent is: 



2pi{^ + ipi)^ [pi + p2)iV^ + ipi){VoJ + ip2) 

{UJ2 ~ OJ) + 



( 



G{lj) = det G{lj) 



{UJ2 - LO) + 



7rA2 



2p2{^/uj + ip2Y 



\ {pi + P2){V^ + ipi){V^ + ip2) 



2p2(v^+«/02)^ / 

^A^ 

(pi + P2)iVi^ + ipi)iVuj + ip2) 

7rA2 



The determinant detG(a;) is 

{detG{u))-^ = 



LOl — bJ 



ttA^ 



2pi{^ + ipiY 
7rA2 



CJ2 ^ + 



2pi{^ + ipiY 



ttA^ 

2p2{y^ + ip2Y 



[Pl + P2){V^ + ipi){V^ + ip2) 



(43) 



(44) 



(45) 



Here we can formulate necessary conditions for the analyticity of the function Gj^}, on the first sheet: 
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These conditions are definitely satisfied in the weak coupling regime, because coi, pi, and A are independent parameters 
and for any fixed toi and pi, in the limit A — > (46) becomes 

1. wipf > 0, i = 1,2 

2. LU1UJ2P1P2 > ■ 

which is obviously true as Uj and pj are positive for any i. 

In order to find out the analytic structure of G'(w), we analyse the poles of the determinant: 



(detG(cj))-i ^ 



(iVl + x^){x + pif - ^ 
2pi 



{l02 + X^)ix + p2f ~ ^ 



ttA^ 



(Pl +P2) 



0, 



(47) 



where we substitute -^/cJ = ix. This is an algebraic equation of 8"^ degree with real coefficients, so all the roots of this 
equation are either real or complex conjugated pairs. All roots are on the second Riemann sheet, and there can be k 
(A; = ... 4) pairs of complex conjugated roots and (8 — 2k) real roots corresponding to virtual states, i.e. negative 
energy states on the second sheet. For weak coupling A — > 0, the third possibility is realized and we have two pairs of 
complex conjugated roots Zj, z*, which can be evaluated perturbatively as: 



2\4 



vr^A 



1 



2pj [ujj + p2)2 ^ ip.y y (^^. _ ujk){pi + P2)2(V^ + ipk) 

+ 0(A6), j = 1,2, k^j. 



For the weak coupling regime the expressions for the real and imaginary parts of Zj are: 

ttA^ - p] 



1] 



Rez," 



-Imz, 



2pi K +P])' 



+ 0(A^), j = 1,2. 



Neglecting the integral term in the representation (|39|), we can write: 

A{t) 



+ 0(A4), j = 1,2. 



fc,fc' = l,2 j = l,2 

k,k' = l,2 
2 2 / 



1 YpAPj 



27rRe(aia2)e 



(pi + P2)(pj - iy/uJJ){pi - iy/ujj)iujj - UJl) 



+ 0(A4), l^j, 



where 



UJl — UJ2 



(48) 



(49) 



We would like to notice that both expressions (^8|) and (49) contain the term \/{LJk — uJi) and, therefore, cannot be 
directly used in the case of degenerate levels in the initial Hamiltonian H^. Also, the case of the continuous spectrum 
of requires a special consideration. 

For the initial conditions oi — 1^ a2 — 0, the survival amplitude ( ^ ) does not have any oscillations. However, such 
oscillations appear in the next order A** in expression (|4^). The survival probability p{t) in the lowest order of A^ can 
be now expressed as: 



p{t) = \A{t)\' = \\a,\'e 



2^-71* 



|a2pe-^^*e-2^.t|2^ 



(50) 



We illustrate the possible behavior of the survival probability in Fig. 2. One can see that depending on the initial 
conditions, the decay can either mimic the behavior of the usual one level model or display considerable oscillations. 
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V. AT-LEVEL MODEL 



In the weak coupling regime we can also analyze the A'^-level model with an arbitrary formfactor Using the 

representation (^), wc find 

N N N 

(detG(c^))-i = J]K-c^)-A2^4fcH n(^™-^)+0(A4), (51) 

k—1 k—1 m^k 



where 



hiiuj) = / duj — — . 

' uj' — uj — lO 







The zeros of this expression give the position of resonances: 

Zk = LUk - X^hkiuJk) + OiX*) ^Uk-ilk, j - 1 ... TV. (52) 
In the first non-trivial order of the perturbation theory with respect to we have: 

Wfc = Wfc, Ik = 7rA^/^(cjfc). 
The partial resolvent G can also be calculated: 

Gkk'i^) = K -UJ- X^hk-y^y^ Skk' + OiX"). (53) 
From this representation we obtain the expression for the residues (Eot) : 



r 



kk' 



Skk'Skj+0{X'). (54) 



We derive the survival amplitude (p9[) in the first non-vanishing term of the perturbation expansion with respect to 
A^: 



A{t) = |afcpe~*'^'=*e-'^^ f"'-'^''^'. (55) 

k=l 

In the case of the iV-level model, the behavior of the survival probability is much more complicated than in two level 
model. In order to illustrate this, we plot in Fig. 3 few examples of the survival probability corresponding to different 
initial conditions for the three level model with different parameters. In this case, the behavior is not necessarily 
"self-similar" even for the very slow decay. One can see that our curves reproduce fairly well the experimental results. 
Hence we can suggest here an explanation of the results which does not refer to the semiclassical description invoked 
in paper Q. Namely, the initial laser impulse creates in the system Nal a wave packet which is a superposition of 
(many) excited states. Then each excited state decays independently while the common survival probability ( [ss] ) 
exhibits a complicated behavior similar to one of Fig. 3 and Figs. 3, 4 in paper ||]. 

In fact, the interference of many decaying states can drastically change the decay patterns. In this case, the decay is 
equally defined by both the parameters of the system (energies, widths) and the distribution of the initial wavepacket 
Uk- Therefore, the decay profile may mimic different non-exponential functions. Such a behavior is illustrated in 
Fig. 4, where we plot the decay of the {2N + l)-level system: 

k 

lk=l = const, Ljfc = cjo + = ■■■N. (56) 

The initial distribution is chosen to be the Gaussian one: 

ak = ^^, afc=exp(-(A)2)^ k^-N...N. (57) 

From Fig. 4 we can see that the initial decay is almost independent of N . The duration of the initial decay is much 
shorter than the decay time Tdec — I/7. The number of repeated peaks decreases as N increases. Already for 10 
levels, a rather small number for molecular systems, the second peak is very far from the region of the initial decay. 
While the values (|56|), ( |57| ) are chosen quite arbitrarily and can only be used for illustrative purposes, we would like 
to notice that the excitation process in experiments similar to is usually well-defined and well-reproduced. Hence 
the initial wavepacket may also be well-correlated. 
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VI. CONCLUSIONS 



Dissociation processes like the dissociation of Nal, which is a kind of tunneUng/decay process, may be described by 
the simple quantum mechanical model of the interaction of the iV-level discrete spectrum with the continuous spec- 
trum. Already the model with two levels displays decaying oscillations of the survival probability in the "exponenial" 
era, while one-level model exhibit the purely exponential decay. The amplitude of the oscillation is determined by the 
initial state, which is a superposition of two excited levels. The model with tree levels may illustrate qualitatively the 
experimental curve of the Nal dissociation. In the iV-level system, the decay is equally defined by both the parameters 
of the system and the distribution of the initial wavepacket. 
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APPENDIX A: TIME EVOLUTION IN TERMS OF GAMOV VECTORS 

By analytic continuation to the second sheet, we obtain the extended distributions G^j(a; ± iO) and — Zfc] + 
defined as functionals, which act on a suitable test function h{Lu) as: 



dLuh{uj)Gti{uj ± iO) = / h{uj)Gki{io±iO) 



(Al) 



OO 

J dujh{Lj)Gkiiio±iO) + 27riY^ J h{uj)Gki{uj ±iQ), 



duj- 



-p uj — ^j 



(A2) 



The contours F and C^^ are presented in Fig. 1. Using ( Al) and ( A2), we obtain from ( p2[ ) the Gamov vectors p8| , p9| , p8| 
in the form 



N 



k,l=l 
N 



k,l = l 



'-31 + 



N I 

A E ML^)GUu + ^0) J du'-^i^^\co') - \k) 



k,l=l 
N 



A E h{Lo)Gki{Lo-iO)\ I Alo' — 

J 1 I ^ 



, A/,(c.') 



kl=\ 



-I- ZO 



(A3) 
(A4) 
(A5) 
(A6) 



We recall that r;[; is the residue of Gki (w -f iO) at the pole Zj . The normalization constants Nk are: 



N 



k,l,m,n—l 



J 

ran 



X\fkiu;)fm{co) 



^km I / 2 

\uj - z., 





(A7) 



The Gamov vectors ( |A3| - A6) are left and right eigenfunctions of the extended Hamiltonian, which can be written as 
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J 
The Gamov vectors form a biorthonormal set: 



which is complete. The completeness follows from the extension of (14): 

N 



k=i i 

The time evolution of the vector \k) in the new extended representation is 

oo 

J 

Using we express the transition amplitude 

N 

J i,i'=i 
N °9 

+ J2 dcue-''^'X^fi{uj)fi,{uj)Gti{u: + iO)Gk'i'{Lo-tO). 
i,i'=i{ 

The integral terms of ( Al2| ) can be rewritten in the form 



1 
27ri 



(A8) 



(A9) 



(AlO) 



(All) 



(A12) 



(A13) 



Taking into account that G^^,, + iO) implies integration along the contour F* , which goes to the second Riemann 
sheet below all the singularities of Gkk'{<jJ + iO), we obtain the transition amplitude in the form 



AT 



k.Lk'.l' = l 



dcje-^'^'Gkk' , 



(A14) 



Ci 



which must coincide with the result obtained using Friedrichs solution (pQl). In order to fulfill this requirement the 
following formula must hold: 



N 



■0)' kl' k'V 



' kk'- 



(A15) 



k,Lk',l' = l 
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Figure captions 

Fig. 1. The contours of integration C and Ci. 

Fig. 2. The survival probability p{t) for the two level model. The parameters are chosen to be 71 = 72 = 10~^, 
u)i = 1.0, LOi = 1.06. The initial conditions are ai = 0.5, 02 = (the solid line), ai = 0.5, 02 = 0.2 (the long-dashed 
line), ai — 0.5, 02 = 0.5 (the short-dotted line). 

Fig. 3. The survival probability p(t) for the three level model. The parameters are chosen to be 71 = 72 = 73 = 10~^, 
oj\ = 1, the initial conditions are ai = 0.3, 02 = 0.5, 02 = 0.3. The energies are UI2 = 1.04, ws = 1.15 (the long-dashed 
line), u!2 = 1.06, u!^ — 1.15 (the solid line), uj2 = 1.064, = 1.15 (the short-dashed line). 

Fig. 4. The survival probability p{t) for the A^-level model. The parameters are chosen to be jk = 10~^, 
Wfc = Wo + kAu)/N, k = —N...N, where wq = 1) = 0.1. The initial conditions are ak = exp(— (fc/Ai")^), 
k = —N...N. The results for AT = 2 (the solid line), N = 3 (the short-dashed line), and N = 5 (the long-dashed 
line) are presented. 
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